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ABSTRACT 


Local anode heat flux in a quasi-steady MPD thruster 
is measured by thermocouples attached to the inside surface 
of a shell anode. Over a range of arc currents J from 5.5 
to 44 kiloamperes and argon propellant mass flows m from 1 
to 48 g/sec , with the ratio J /m held constant, the fraction 
of arc power deposited in the anode is found to decrease 
with increasing arc power. Specifically, this anode power 
fraction decreases from 50% at 200 kW arc power, to 10% at 
20 MW. In an effort to account for this functional behavior, 
the current density, plasma potential, and electron tempera- 
ture in the plasma adjacent to the anode are measured with 
probes, and the results are used in a theoretical anode heat 
flux model. The model asserts that energy exchange between 
electrons and heavy particles in the plasma near the anode 
occurs over distances greater than the anode sheath thick- 
ness. Consequently, a contribution from interelectrode 
potential augments anode fall voltage, electron temperature 
and work function contributions to the anode heat flux. 
Calculations of anode heat flux using the model agree with 
the direct measurements, and the decrease in anode power 
fraction at higher arc powers can be attributed to the 
smaller mean free paths in the interelectrode plasma. 
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CHAPTER 1 

THE ANODE POWER QUESTION 


1.1 Introduction 

Electromagnetic thrusters accelerate ionized propellant 

by the interaction of an electric current driven through it 

with a magnetic field. This class of electric propulsion 

promises a combination of high specific impulse (1000-5000 

4 2 

seconds) and thrust density (over 10 N/ra ) demanded by many 

( 2 ) ... 

deep space missions that is unattainable with chemical 
propulsion, ^ electrothermal ^ ® ^ or electrostatic 
thrus ters . ^ ^ ^ ^ In electromagnetic thrusters using cur- 

3 

rent levels below 10 amp, an external magnetic field must 
be supplied to produce significant thrust; at higher 

currents, the self-induced magnetic field is of sufficient 
magnitude that the external field may be avoided. it 

is with this latter, self-field class of electromagnetic 
thruster that this study is concerned. 

The particular type of electromagnetic thruster examined 
here is an outgrowth of earlier attempts to extend the opera- 
ting range of thermal arc jets by increasing their power dens- 
ity. A typical thermal arcjet with a specific impulse of 
about 1000 seconds and a thrust of about 1 Newton has an arc 
power of about 30 kw. Attempts to improve specific impulse 

and thrust by increasing input power caused severe heat load- 

( 9 ) 

ing of the anode. 

Although the performance capabilities of the arcjet had 
been improved by the use of magnetic fields, the electrodes, 
especially the anode, were still being heated excessively. ^ ^ 
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However, working with a conventional thermal arc jet with a 

central cathode and a coaxial anode, Ducati showed that 

drastically reducing the mass flow caused the establishment 

, (17) 

of a diffuse axisymmetnc discharge. In this mode, the 

arc current could be increased and the power range in the 
steady arcs enlarged by an order of magnitude. With arc 
currents in the kiloampere range, the induced magnetic fields 
were sufficiently large that external magnets were unnecessary 
This thruster, the self-field MPD arc, shows the potential 
required for .many advanced deep space missions. 

A schematic of such an MPD arc is shown in Fig. 1.1. 

The current, which flows from anode to cathode, has axial 
and radial components and produces an azimuthal magnetic 
field. The arc current passes through the injected propel- 
lant stream, ionizes the propellant and interacts with the 
induced magnetic field to accelerate propellant mass to pro- 
vide thrust. The interaction of the radial component of arc 
current with the azimuthal magnetic field produces an axial 
acceleration ("blowing" - ) , while interaction of the axial 
component of arc current with the magnetic field compresses 
the propellant to provide a reaction force on the cathode 


surface (“pumping"). In addition to these electromagnetic 
components of thrust, there may be some electrothermal con- 
tribution recovered in the expansion of the high-enthalpy 
propellant heated ohmically in the arc chamber. 


The electromagnetic thrust in the MPD arc varies as the 

(20 21 ) 

square of the arc current. # ' This suggests that the arc 

be operated at the highest possible current levels. However, 
steady high-current operation is limited by tolerable elec- 
trode heat transfer as well as by power supply capabilities. 

To relieve this difficulty, the arc may be operated in a pulsed 

(2223) 

"quasi-steady" mode. # In quasi-steady operation, electri- 
cal power is supplied to the thruster in pulses which are long 
enough that steady processes prevail during the discharge 
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period, but short enough that the average energy deposited 
in the electrodes is tolerable. The quasi-steady operation 
also facilitates laboratory study of the MPD accelerator. 

It allows detailed study throughout the plasma using un- 
cooled diagnostic probes. 

The quasi-steady thruster used in this investigation 

has been extensively studied at Princeton. A discussion 

(22 23) 

of its development is discussed by Jahn and Clark, ' 

( 22 ) ( 20 ) 

the performance described by Clark, Cory and 

.Bruckner, (24) an ^ .the .mechanisms involved investigated by 

( 25 ) . . 

Boyle and Jahn. Examination of the cathode region was 

discussed by Turchi^^) and anode phenomena were first 

(27) 

examined by Oberth. The accelerator was designed to 

study the physics of MPD acceleration, but results apply to 
the development of the MPD arc jet as an electromagnetic 
thruster . 


1 . 2 Anode Power 

The efficiency of electromagnetic thruster operation at 

powers below one megawatt is dominated by anode losses, which 

( 28—30 ) 

may exceed 50% of the total arc power. ' Because of 

this and the electrode heat transfer problem, considerable 
attention has been paid to the evaluation and understanding 
of anode power in these devices. 

The anode power in steady thrusters is usually assessed 

from the temperature rise of a known mass flow rate of anode 

(31) 

coolant. In a variant of the technique, the anode power 

of a radiatively cooled 50 kW external field accelerator was 
estimated to be 20-40% of the arc power from the heat ab- 
sorbed by the magnet housing coolant . ^ ^ Another investiga- 
tion used the temperature rise of the coolant of a segmented 

anode to estimate local anode heat flux in a sub-megawatt 

(32) 


MPD arc. 
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Existing theories of anode heating are based on energy 
balances over the anode sheath in a plasma free of magnetic 
field and in thermal equilibrium. (33,34) The results have 
been verified in high-pressure arcs operating with currents 
from 1 to 500 Amps and current densities from 1 to 500 A/cm 
with arc powers to 50 kW. 31,35,36) However, these theories, 
as well as those dealing with the flow of an initially un- 
ionized gas presumed to ionize over the anode/^'38, 39,40,41) 
do not strictly apply to relatively low pressure discharges. 
Nevertheless they have been used to infer the anode power in 
the MPD arc . 

In earlier experiments in this laboratory, it was shown, 

from measurements of local plasma properties, that the anode 

power fraction in an MPD thruster decreases as arc power is 
(27) 

increased . This result, based on indirect measurements, 

suggests that a fractionally larger portion of the input power 
is available for the useful production of thrust. Consequent- 
ly, detailed study of anode heat flux and anode power is 
warranted. 

In the investigation reported here, local anode heat 
■^1^^ and power in the MPD arc are measured directly over an 
extended range of arc powers up to multi-megawatt levels, 
and a theory of anode power which is applicable to the MPD 
arc is introduced. The investigation is carried out on an 
apparatus described in Chapter 2, and the direct measurement 
of local anode heat flux is shown in Chapter 3. A study of 
the anode region of the discharge is presented in Chapter 4 
and the mechanisms involved in anode heat flux are detailed. 
Chapter 4 also includes an explanation for a measured de- 
crease of anode power fraction with increasing arc power ob- 
served in Chapter 3. The implications of the work to high 
current discharge devices, and the MPD arc in particular, are 
discussed in Chapter 5. The appendices which follow detail 
several of the arguments and calculations used. 



CHAPTER 2 

THE MPD ARC JET APPARATUS 


2.1 Introduction 

An MPD arcjet apparatus, designed and built specifi- 
cally for this study, is described below. Included are 
descriptions of its discharge chamber, anode, propellant 
injection system, power supply, gas-triggered switch and 
vacuum facility. The operation of the device and its 
voltage— current-mass flow characteristics are also in- 
cluded c 

2.2 The Discharge Chamber 

The apparatus is shown schematically in Fig. 2.1 = The 
discharge chamber is a 12.7 cm diameter by 5 cm deep cylin- 
drical cavity formed by a Plexiglas insulating ring 0.64 cm 
thick. It is backed by a 1 cm thick Plexiglas plate sur- 
rounding a central, conical, 2%-thoriated tungsten cathode 
2.7 cm high with a 1.9 cm base diameter. Around the cathode, 
arranged at the corners of a regular hexagon, 2.5 cm from the 
axis, are six propellant injectors. They are made of nylon, 
have 4 mm orifices and extend 6 mm from the backplate surface. 
The coaxial anode is 1 cm wide with a 10*2 cm diameter ori- 
fice. 


2.3 The Shell Anode 

The heart of the MPD arcjet under study is the 1 cm wide 

aluminum shell anode of 1 mm wall thickness. The 

anode is masked with mylar insulating tape to expose the 

2 

same surface area for current conduction (^260 cm ) .as the 
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(27) 

usual solid anodes employed in thus laboratory. The 

1 mm anode wall thickness is a compromise between ade- 
quate strength and the desire to use a 1-dimensional . 
analysis of heat flow through the shell for the anode 
heat flux measurements. The anode is set at ground 
potential by connecting it to a grounded pipe with a 
#16 ga (1.3 mm diam) insulated copper wire. 


2.4 Propellant Metering and Injection 

Argon propellant from a gas bottle supply is fed 

to a 300 ml cylindrical reservoir, 5 cm o.d. x 24 cm 

overall length. The downstream end of the cylinder is 

fitted with an electrically actuated valve which leads 

to the injectors. The upstream end is connected by one 

side of a "T" to the supply via a metering valve, while 

the other side of the "T" connects the cylinder, through a 

valve to a vacuum pump. For operation with plenum pressures 

above atmospheric upsteam of the valve (for argon mass 

flows greater than 2 g/sec) adjustment of bottle pressure meters 

the propellant , and for operation requiring subatmospheric 

pressures, the cylinder is isolated from the gas supply 

bottle and its pressure is adjusted by the valve and 

vacuum pump combination. The reservoir is large enough 

to maintain essentially constant mass flow of propellant 

throughout the complete discharge. The valve is opened 

in about 2 msec using an electronic trigger (Fig. 2.2),^^ 

(Fig. 2.2),^ 42 ^ by suddenly overdriving the 12 volt rating 

of the valve solenoid by a factor of 10. Constant mass flow 

is established after the valve response risetime of less 
(43) 

than 3 msec. 
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Mass flow is calibrated against reservoir pressure 
using a Varian "Milliter" (T,M.) fast ionization gauge 
mounted on axis 1 cm from tbe cathode tip, to determine 
the duration of the constant mass flow and a McLeod gauge to 
determine the pressure rise in the vacuum tank, for a cold 
flow (no arc discharge). The experimentally derived pro- 
portionality between pressure and mass flow is constant, 
demonstrating that the flow is choked at a constant 
area. ^ 100 ^ The accuracy in the calibration is t5 %, 
determined by the deviation of the calibration data 
from a straight line. 

The mass flow into the thruster is limited on the 
low extreme (1 g/sec) by proper seating of the valve and on 
the upper (50 g/sec) by the gas reservoir pressure which 
holds the valve shut. 

2.5 Power Supply 

A rectangular, nonreversing, current pulse is sup- 
plied to the discharge chamber by a 160 kilojoule, 3300 \xF , ■ 
capacitor bank with a total stored charge of 24 kA-msec, 
i.e. Coulombs, when charged to 8 kV. . * 44 ^ This power supply, 
an L-C ladder network (Fig. 2.3) consists of 119 series 
inductors and 120 shunt capacitors (G.E. Capacitor Unit, 

Cat. No. 14F641) . The capacitors are housed in a plywood rack 
30 units across by 4 units high. Since each capacitor 
has its own compartment in the rack, adjacent units are 
protected in case of failure^ One terminal of each capaci- 
tor is connected to the inductance line using a 10 cm long 
fuse of #16 ga (1.3 mm o.d.) tinned copper wire enclosed 
in a nylon sleeve. The other terminal, and for safety the 
metal capacitor case, are connected to ground. The inductors 
fabricated from 3.8 cm wide x 0.32 cm thick copper strip, have 
three 10 cm i.d. x 6.4 cm pitch turns between capacitors. 

The end of each tier is connected to the next by an 
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inductor designed to maintain the flatness of the current 
pulse. 

The pulse is made nonreversing by matching the source 
impedance (ladder-network) to the load impedance (MPD arc) 
using a ballast resistor for each tier in the form of an 

electrolytic cell consisting of two copper plates of 2000 

2 . . 
cm surface area on each side immersed in a 60 liter vat 

filled with 0.58 molar cupric sulphate solution. 

The tiers can be connected in different arrangements. 
For example, in series for a- 900 psec pulse (maximum cur- 
rent: 26 kA ) 7 a pair of two series-connected tiers in 
parallel for a 450 psec pulse (maximum current: 52 kA) ? 
or all four tiers in parallel for a ,225 psec pulse (maxi- 
mum current; 104 kA) . 

2.6 The Gas-Triggered Switch 

The power supply is connected to the accelerator 
through a gas- triggered switch shown schematically in Fig. 
2.4. The switch is formed from two 13.8 cm diameter by 
1.9 cm thick stainless steel plates separated by a 12.7 cm 
i.d. x 5 cm long Plexiglas ring with a 0.64 cm wall. The 
assembly is housed in a cylindrical aluminum case, which 
also forms part of the electrical return path to the power 
supply. The switch is prepared for a pulse by being evacu- 
ated to about 10 mt^rr absolute pressure via a hole in one 
plate, through a vacuum dump tank, by a mechanical pump. 
Nitrogen gas, passed into the switch through an electrically 
actuated valve and a diffuser mounted in the opposite plate 
triggers electrical breakdown allowing current to flow. 

2.7 Vacuum Facility 

The vacuum facility is shown in Fig. 2.5. The accel- 
erator is mounted in the end of a 45 cm diam x 76 cm 

















FIGURE 2.5 
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diameter access ports, one on tine end and the other on the 

-2 . 

side* The background pressure of about 10 mtorr as estab- 
lished by using a 5 cm oil diffusion pump with a 425 liter/ 
minute (free air) mechanical forepump* A 2*5 cm thick Plexi- 
glas slab sealing plate with 0.64 cm thick aluminum backing for 
support is fitted with a manual probe mount and traverse 
mechanism capable of moving diagnostic probes radially and 
axially throughout the discharge without opening the system* 

For safety, the entire assembly is mounted on a table and 

housed- in - a - "Plexiglas" case - ," behind a Plexiglas wall* 

2.8 Operation of the MPD Apparatus 

To operate the thruster, the capacitor bank configuration 

is arranged and the power supply connected to the gas switch* 

-2 

The vacuum tank backpressure is established, usually at 10 
mtorr. The capacitors are then charged to the voltage appro- 
priate for the current desired. Closing the firing switch 
triggers the propellant valve and initiates an electronically 
delayed pulse which trips the gas switch valve* Precise delay 

of the discharge is necessary to achieve proper arc opera- 
(23 ) 

tion. x ; If the delay is too long, the large amount of pro- 
pellant mass injected prior to the discharge increases the 
backpressure in the vacuum tank and compromises the simulation 
of a space environment. If the delay is too short, the dis- 
charge is initiated before the mass flow is steady, resulting 
in time-varying discharge processes. After the completion of 
the discharge cycle, the system is evacuated in preparation 
for another pulse. 

2.9 Current-Voltage-Mass Flow Characteristics 

The apparatus is operated with arc currents from 5*5 to 
44 kA and injected argon propellant mass flows from 1 to 
48 g/sec a 
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The current is measured using the passively integrated, (46) 

output of a 340-turn, calibrated Rogowski loop encircling the 

(47) 

discharge chamber. 1 The arc voltage is simultaneously 

3 8 

monitored with a Tektronix 6013A, 10 times attenuating, 10 ft 
impedance, voltage probe, mounted at a cathode connection. 

A typical response, for 900 psec operation at 16 kA arc cur- 
rent and 24 g/sec injected argon mass flow is shown in Fig Q 
2.6. The current rises in about 10 psec while the voltage 
has an initial transient of about 150 psec after which it 
holds a steady value of about 80V, until the discharge ceases. 
The term "quasi-steady" has been applied to this later period 
during which relatively constant operation ensues. 

Terminal voltages of the discharge for various arc cur- 
rents and argon mass flows are shown in Figs. 2.7 and 2.8, v 
respectively. None of the characteristics is a straight line: 
between two apparently curved sections there are jumps notice- 
able in all the characteristics. The terminal voltages agree 

with those of a like-geometry thruster mounted in a 1 m diam 

( 22 ) 

x 2 m long Plexiglas vacuum tank. ' 
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CHAPTER 3 

ANODE HEAT FLUX AND POWER 


3 o 1 Introduction 

The direct measurement of anode heat flux and power 
are described in this chapter. The calorimetric technique 
employed uses thermocouples attached to the inside surface 
of a hollow shell anode in a series of experiments over a 
range of arc powers from 200kW to 20MW, using arc currents 
from 5.5 to 44 kA with argon propellant mass flows from 
1 to 48 g/sec. The reader is referred to Appendix E for 
a compendium of all the results. 

3c2 Anode Temperature Rise 

The MPD thruster anode has a 1 cm wide anode with a 

10 cm diam orifice. These dimensions permit use of a shell 

anode with a 1 mm wall (Fig. 3.1). The anode is internally 

fitted with twelve #34 ga (0.16 mm diameter) copper-constantan 
( 48) 

thermocouples. The thermocouples are welded to the in- 

side surface within a 15° arc, usually radially opposite a 
propellant injector. Two are in the chamber side (A, B) , 
five in the lip region (C,D,E.F.G) and five in the downstream 
face c These thermocouple junctions, formed by arc welding 
wire ends directly to the inside surface of the shell, are 
the "hot” junctions for the measurement 0 

To measure changes in their potential, every thermo- 
couple inside the anode is supplied with its own reference 
'•cold" junction thermocouple outside the anode. The cold 
junctions are individually made by welding the tips of 
thermocouple wires together and then electrically insulating 
the resulting bead with epoxy 0 All twelve cold junctions 


! 
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are bundled in a Pyrex tube sealed at one end and the 
package is immersed in a half-liter Dewar flask filled 
with ice water to maintain a constant temperature. 

Each thermocouple pair produces an output voltage 
which changes with hot junction temperature. This 
voltage is monitored using the detection circuit seen in 
Fig . 3.2. The circuit incorporates a 10 kn, 22 turn bias 
resistor to null ambient output due to room temperature. 

The circuit differentially amplifies by lOOOx the net 
voltage change, which in some cases is only tens of micro- 
volts . The input impedance of the preamplifier is 1000ft , 
compared to the 25ft output impedance of each thermocouple 
pair. 

The amplified output signals of all 12 thermocouple 
pairs appear simultaneously on Tektronix 555 Dual-Beam 
oscilloscopes equipped with multiple beam Type M and 
1A4 plug-in pre-amplifiers and are recorded on Polaroid 
film. The amplifier for each thermocouple pair is cali- 
brated over the frequency range from 10 to 1kHz using a 
sine wave generator and the entire anode-thermocouple- 
amplifier-oscilloscope calibration is checked by immersion 
of the anode in a fixed temperature water bath. The time 
response of the system is 3 msec. It is measured by plac- 
ing a cube of water ice on one of the hot junction thermo- 
couples and recording the time required for the output 
signal to drop from room temperature to a constant level 
for the temperature of the ice. 

Connections between dissimilar conductors anywhere in 
the circuit produce potential differences changing with 
temperature. Therefore, the introduction of an error re- 
sulting from differentials in ambient temperature between 
the connections of the hot and cold thermocouples in each 
pair is precluded as shown in Fig. 3.2. The constantan 
leads of the hot and cold junctions are welded together and 
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the copper wires from each thermocouple pair are soldered 
to the copper input buss of the amplifier. 

The MPD discharge pulse generated substantial electro- 
magnetic noise. To eliminate it, inductance of the circuit 
is minimized by using sub-miniature components in construc- 
tion and by twisting all leads in pairs; and capacitive 
coupling to the high voltage discharge breakdown potential 
is reduced by covering the detection circuit and all of its 
connections with a copper shield grounded to the anode. 

Response of the detection system to a 900 p sec long, dis- 
charge with 16 kA arc current and 6 g/sec injected argon 
propellant is illustrated in Fig. 3.3. Thermocouple signals 
do not begin until the initiation of the discharge, shown 
on the 5 msec/div time scale G When the discharge starts, 
electromagnetic noise distorts the baseline and the output 
of the thermocouple detection circuit. However, the noise 
subsides during the 10 msec required for the temperature 
increase of the unheated side of the anode shell to be 
nearly complete, when the output signal becomes almost con- 
stant. Subsequently, the outputs of cooler thermocouples 
(A,B,C,H,K,L,M) continue to rise, while those of relatively 
warm thermocouples begin to droop. Data are examined 20 msec 
after initiation of the discharge, when the noise has sub- 
sided and heat conduction along the shell is still negligible 
(Appendix A) . One of the thermocouples, P, is situated under 
mylar electrical insulation at the periphery of the plasma- 
exposed region of the anode. A temperature rise at its loca- 
tion is due only to current conduction along the anode and its 
response is null for the first 10 milliseconds following the 
discharge c Therefore, any temperature rise due to Ohmic dis- 
sipation in the aluminum anode is neglected. The signal at 
later times is due to heat conduction from the adjacent anode 
region which is exposed to plasma. 

Operation is repeated 3 to 10 times at each current 
and mass flow condition, the thermocouple output voltage 
is averaged, and the temperature rise is then calculated 
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from the thermocouple calibration. The temperature rise 
due to a 900 psec long discharge at 16 kA arc current and 
24 g/sec mass flow is shown in Fig. 3.4, The abscissa 
represents a linearly rolled-out anode: it traces a line 
of constant azimuth along the conduction surface . The 
axis runs from the chamber insulator (shown as a cross- 
hatch) , around the lip (with midplane at 0.0 cm) almost 
to the boundary of the conduction region on the anode face 
(at about 4.5 cm) . The temperature rise peaks at about 
20°C at the lip and drops to less than 0.5°C on the down- 
stream face and to 0 . 1°C on the chamber side. The error 

(49) 

bars shown represent one standard deviation of the sample. 

The temperature profiles measured are all peaked markedly 
at the lip. The values there are more than an order of 
magnitude higher than those on the face or inside the dis- 
charge chamber. 

3.3 Local Anode Heat Flux 

The local anode heat flux during the quasi-steady 
discharge is derived from the measured temperature rise on 
the inside of the anode shell. Heat supplied to the 
plasma-exposed side of the anode is presumed to be con- 
ducted through the shell wall, increasing the temperature 
of the inside surface. Once the temperature rise on the 
inside surface is measured, and the heat flux problem for- 
mulated, the local anode heat flux can be derived. 

The assumptions underlying the analysis follow, 
a) Heat flux to the anode is constant during the discharge 
and zero at other times. Initiation and termination 
transients are negligible for quasi-steady discharge pulses 
larger than 450 psec. Thus, all heating is steady and 
confined to the discharge period. Indeed, a more detailed 
analysis is not warranted since fluctuations in the heat 



“TJ 


> 

U 
ro 

m 

o 

o 


O 

c 


OJ 





ANODE TEMPERATURE PROFILE 
16 kA , 24g/sec 



flux during tlie discharge are smoothed out in the con- 
duction process through the anode shell* 

b) Heat conduction along the shell can be neglected in 
comparison with conduction through the shell, since tem- 
perature gradients along the shell produce negligible 
heat fluxes compared to those through the shell (Ap- 
pendix A) . 

c) The anode is assumed to be locally flat everywhere, 
since the anode is flat on the face and discharge chamber 
sides and the radius of curvature of the lip is much 
greater than the shell thickness. 

d 0 Ohmic heating of the shell due to current conducted 
along it can be neglected D This is verified by the 
negligible response of thermocouple P located under the 
mylar anode masking, immediately following the discharge 
and by analysis of current conduction through aluminum 
slabs . 

e) The temperature of the anode surface does not become 
large enough for radiated heat flux to be significant. 

f) The anode does not ablate and therefore power to erode 
the anode is neglected. Local spot attachment of the arc 
is also neglectedc 


The heat conduction analysis uses a time-averaged 
local anode heat flux defined at each thermocouple 



o 


where ) is the actual anode heat flux, which may be 

a function of time? T is the discharge time? and t f 
usually equal to the discharge time, is the time at which 
all heat flux to the anode ceases. This formulation sup- 
presses any time-varying details of the heat flux during 
the discharge. 
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If an annular anode section is associated with each 
thermocouple, then heat conduction through the shell wall 
can be analyzed as a series of one-dimensional heat flows 
through those slabs. For each section, the temperature 
rise on the unheated side of the shell is known, and the 
heat flux is assumed constant for its duration, therefore 
the heat flux problem can be solved in reverse: What 

heat flux produced the temperature rise? 

The result can be determined from the following re- 
1 at ion— ( Appendix -A-)- 




( 3 0 2a) 


(3,2b) 


for the temperature rise, 0 (t) on the unheated side of 
a uniform slab of thickness Si, at time t>T (the pulse 
duration) due to a heat flux q uniform over the surface 
of the opposite side begun at t = 0, For aluminum, ^ 
c = specific heat = 0.903 joules/ (g-°C) 

^ = density = 2.70 g/cm^ 

k^ = thermal conductivity = 2.01 W/cm-°C 

ae = v<? c = thermal diffusivity = 0.823 cm /sec 

The computed temperature rise per unit heat flux 
^ (t) for the actual anode wall thickness of 1.03 to 1.12 
mm and a 900 psec pulse duration is shown in Fig. 3.5a* 

The rise begins at zero and increases within 10 msec to 
about 3 C/ (kW-cm ) . This value then remains constant 
since, by assumption, the slab retains all the heat sup- 
plied. The appearance of this analytical result agrees 
with the experimental example shown in Fig. 3.5b. There- 
fore, local anode heat flux is computed by dividing the 
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experimentally measured temperature rise due to a dis- 
charge pulse by <^>(t) from Eq. (3.2) . 

The anode heat flux during MPD arc operation with 
16 kA total current and 24 g/sec argon propellant mass 
flow, computed as described above, appears in Fig. 3.6. 

Due to the small variations in shell thickness along the 
anode (1.02 to 1.13 mm), there is a slight difference 
in the shapes of the temperature and heat flux profiles 
on the surface. However, just as the temperature pro- 
file is peaked at the lip, so is the heat flux. The 

2 

maximum heat flux is 6000 W/cm over the lip region just 

downstream of the anode midplane. The heat flux drops 

,2 ,2 
to about 100 W/cm on the face and 40 W/cm in the dis- 
charge chamber . 

The magnitude of the local heat flux changes with 
current and argon mass flow, discharge duration, vacuum 
tank back-pressure, and in some cases is different at 
different azimuthal locations relative to propellant in- 
jectors. 

The effect on anode heat flux of a threefold increase 
in current from 8 to 22 kA, while mass flow is held con- 
stant at 24 g/sec, is illustrated in Fig. 3.7. The change 
of current is accompanied by over an order of magnitude 
increase in local anode heat flux. When the current is 

raised, the heat flux maximum in the lip region changes 
2 2 

from 1000 W/cm to nearly 20,000 W/cm and on the face 
the change is from 20 to nearly 400 W/cm . it is in- 
teresting to note that, for a mass flow of 24 g/sec, the 
heat flux on the anode face when arc current is 22 kA is 
about half the maximum heat flux at 8 kA. In comparison 
to operation with different currents, if current is con- 
stant at 16 kA, while mass flow is changed from 6 to 
48 g/sec, changes in the heat flux magnitude are not very 
great. However, the pattern of the change is a more 
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complicated. As shown in Fig. 3.8, the heat flux at the 

2 

lip rises from about 6000 W/cm when the argon mass flow 

2 

is 6 g/sec to a maximum of nearly 20,000 W/cm at 12 g/sec, 

and as mass flow is further increased, the heat flux at 

the lip falls once again. This behavior at the lip is 

accompanied by a monotonically decreasing heat flux over 

2 

the discharge chamber side from about 100 to 30 W/cm 

2 

and over the anode face from 300 to 100 W/cm « 

These changes in mass flow and current may be accom- 
panied by changes in the character of the arc operation. 

For exampTeV arrt'nc r ease' o f — cu r r e nt~at”f ix ed ma ss f l ow 
may be accompanied by the onset of ablation of insulator 
and anode material. The result caused by the ablated con- 
taminants is increased number density, an important 
factor in the explanation of anode power given in Chapter 
4. 

If the pulse duration is 450 psec or below, the 
initial and final transients constitute too large a por- 
tion of the discharge period to be neglected. This is 
seen when the anode heat flux is measured for operation 
at 16 kA x 6 g/sec for different discharge periods. 

As seen in Fig. 3.9, the heat flux maximum shifts down- 
stream and the heat flux to the anode face increases when 
pulse length is increased from 900 to 225 psec. The change is 
small, but demonstrates that as discharge time increases 
the anode heat flux distributes on the anode. Indeed, for 
some cases where the initial and final transients extend 
for the full duration of the discharge, the arc is not 
quasi-steady and increased anode heating upstream is sug- 
gested. Transient effects were precluded by studying quasi- 
steady discharges of over 450 psec duration. 
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A changing ambient backpressure can also influence 
anode heat flux. While other investigations on anode 
power have used background environments with pressures 
as high as 1 atmosphere and only as low as 10 torr^ 31 ' 35 ^, 
this study has concentrated on anode heat flux for an MPD 
arc jet exhausting into a 10 2 mtorr background pressure. 

A comparison is outlined for the arc operated at .16 kA 
arc current and 6 g/sec injected argon mass flow (Fig. 

3.10) u The backpressure is increased from the usual less 

— 2 

than 10 mtorr to 500 mtorr by prefilling the vacuum tank 
w ibh argon gas. On the anode face the heat flux decreases 
only slightly with increasing backpressure, in the dis- 
charge chamber it does not change, and the maximum heat 
flux at the lip decreases by half from 8000 to 4000 W/cm 2 
as the backpressure goes up. These changes are not ob- 
served at the same current when the mass flow is higher. 
Hence, one may speculate that the major effect of the 
higher backpressure is an increased number density in the 
anode region. The results point out that heat flux measure 
ments must be performed in a vacuum environment closely 
approximating space conditions if performance is to mirror 
the results to be expected in space. 

Figures 3.6 to 3.10 have presented the heat flux at 
anode locations along an azimuth opposite gas injectors. 
However, by rotating the gas injectors 30° around the 
discharge chamber axis, the heat flux at locations be- 
tween gas injectors is found. The results are shown in 
Fig* 3.11, where the heat flux opposite injectors is 
compared to the heat flux between injectors for a cur- 
rent of 16 kA with 32 g/sec injected mass flow. At the 
anode midplane (0.0 cm) and upstream into the discharge 
chamber, the anode heat flux between injectors is higher 
than opposite injectors. Downstream it is lower than for 
locations opposite injectors. The azimuthal variation 
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which appears is mild, with the heat flux values differ- 
ing by less than 10% locally. The variation may be a 
result of nonuniformities seen in optical studies of 

the discharge, where the particle number density distri- 

f 241 

bution appears to depend on gas injection geometry 
and the variation decreases with decreased mass flow at 
the same 16 kA current. In any case, for the purpose of 
the measurement of anode power in the MPD arc, once the 
heat flux is integrated over the entire anode surface, 
the influence of the azimuthal variation is small, falls 
within the experimental error bar, and is entirely neg- 
ligible. 

3.4 Anode Power from Direct Measurements 

Once the anode heat flux is determined, it is inte- 
grated over the entire anode surface to yield the anode 

power . The local anode heat flux q is assumed con— 

8 . 

stant over each annular anode surface area A associated 

m 

with each thermocouple: the area extends from points mid- 
way to adjacent thermocouples, except for thermocouples 
closest to insulating boundaries, whose areas extend to 
those boundaries. Therefore, the anode power is 

~ rr> (3.3) 

nr» s \ 

The calculation is insensitive to the exact selection of 
an area associated with each thermocouple, a consequence 
of the existence of large heat flux gradients only where 
the areas are small, i.e., at the lip. Furthermore, if 
the data are examined at any time from 15 to 400 msec fol- 
lowing the 225 to 900 usee discharges, the value of the 
anode power does not change. This verifies that selection 
of 20 msec as a time to examine thermocouple data allows 
the local heat flux to be calculated, while the integral 
of the flux over the anode accurately represents anode 


power . 
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In the following, anode power is reported as a frac- 
tion of total arc power, P^/P^, to emphasize the relative 
importance of this component of the power allocation. For 
example, for the MPD arc operating with a total current 
of 16 kA and an injected argon mass flow of 24 g/sec, the 
anode power of 260 kW represents 20% of the 1,3 MW arc 
power 0 This implies that 80% of the power is allocated to 
other parts of the discharge. 

If the current is held constant, anode power fraction 
peaks with mass flow. Figure 3.12 shows anode power fraction 
maxima occur as mass flow extends from 1 to 48 g/sec, while 
the current is fixed at 11, 16 or 22 kA. When current is 
16 kA, for example, the anode power fraction for a 3 g/sec 
argon mass flow is Chi. The fraction rises to a peak of 
about Ch24 to 16 g/sec and then drops again to nearly 0.1 
as mass flow reaches a 48 g/sec level. 

The peaking of anode power fraction occurs for all the 
currents studied, shifting to higher mass flow as the cur- 
rent level is raised c The search for a correlation factor 

n m 

has led to the consideration of J^/m, the ratio of the 
square of the current to the argon mass flow. This para- 
meter has appeared in various other contexts in MPD arc 
studies in this and other laboratories. It has served to 
correlate arc geometries as well as phenomena such as the 
appearance of fluctuations in arc voltage ^ 51 ^, and the at- 
tainment of a maximum exhaust velocity with increasing 
current for a given propellant mass flow c ^ 52 ^ When the 
anode power fraction P A /P T is related to J 2 /m (Fig. 3*13), 
the peaks of anode power fraction for all currents occur 
j /m between 15 and 20 kA 2 — sec/g„ For values of J 2 /m 
below this, arc operation is quasi-steady and relatively 
free of ablation; if J 2 /m exceeds 16 kA 2 -sec/g, the ablation 
can no longer be neglected. When J 2 /m is larger than 
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16 kA -sec/g the argon mass flow does not represent the 
actual flow rate of particles through the accelerator. 

The particle number density increases as impurities ap- 
pear as a result of the ablation of the Plexiglas insu- 
lator. Thus, the appearance of low anode power frac- 

tion at high J /m may be artificial. 

Since arc voltage depends on propellant mass flow, 
the same arc power can be achieved using a variety of arc 
currents. However, a single operating condition that 
corresponds to a value of arc power is de.fi ned_i.f-j— /m 
is specified. On this basis, the dependence of anode power 
fraction on arc power can be represented uniquely c This 
dependence is illustrated in Fig„ 3.14, where the anode 
power fraction P^/P,^ is shown over a range of arc power p 
from 255 kW, the level of the high-power steady arcjet, 
to 16.5 MW, nearly two orders of magnitude. The data dis- 
played cover a range of arc operating conditions from 
2 ,* 2 

j /m — 5 and 10 kA -sec/g where little or no ablation of 

2 • 9 

chamber insulators occurs, to J/m = 20, 40 and 80 kA -sec/g 
where the increasing ablation cannot be neglected. Illus- 
trated is a decrease of anode power fraction as arc power 
• 2 

increases, even at the highest J /m conditions. The results 

show the trend to decreasing anode power fraction whether 

propellant is supplied independently or, as in vacuum or 
(53 54 } 

ablation arcs, ' the propellant is supplied by the 
ablation of a solid insulating material. 

Figure 3.14 is a summary of the directly measured anode 
power fraction sought in this study 0 There is an overall 
drop in anode power fraction from 50% of the arc power at 
200 kW to less than 10% as arc power is increased to nearly 
20 MW. The important implication is that more power is 
available to the plasma, some portion of which will appear 
as thrust power, yielding more efficient thruster operation,. 
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CHAPTER 4 

ANODE HEAT FLUX MECHANISM 


4.1 Introduction 

One model commonly employed to describe anode heat flux 
in an arc discharge free of magnetic field presumes an elec- 
tron energy balance over the anode sheath. In this 

model, the local anode heat flux, q , is related to the flow 

a 

of electrons of temperature T with a current density j 

e a 

through the anode fall voltage, V , to an absorbing anode of 

a 

work function^; 

V" * ‘V™ 1 " (4.D 

where k is Boltzmann's constant, e is the charge on an 
electron and <3 conv an ^ c 3 ra ^ n are the convective and radiative 
contributions. This model has been successfully applied to 
atmospheric and high pressure arcs where the anode is immersed 
in a dense plasma . ( 3 1 ' 35 ^ However, application of the model 
in this form to the MPD discharges under study here is question- 
able, in view of the much lower number densities and significant 
magnetic induction field which prevails. The effective energy 
exchange mean free paths in the MPD discharge, which are in- 
fluenced by the magnetic field, are much larger than the anode 
sheath thickness. Consequently, while the electron tempera- 
ture term is still valid and the work function still applies, 
the anode fall voltage no longer represents potential energy 
the electrons gain from the electric field and deliver to the 
anode as heat. 

This chapter describes the determination of the properties 
of the plasma adjacent to the anode, and a modification of the 



48 


anode heat flux theory to make it applicable to the MPD thruster. 
Contributions from each term in the model are then evaluated and 
the anode power is calculated. The results agree with the mea- 
surements reported in Chapter 3. Based on this agreement, the 
MPD anode heat flux model is used to explain the measured de- 
crease of anode power fraction with increasing arc power. 

4.2 Anode Current 

The -looa-1- a-node current density is evaluated from experi- 
mental maps of the magnetic induction field in the discharge. 

The magnetic field is measured with a pair of induction coil 
probes. Each probe has a coil of #38 ga (0,01 mm diam) 
enamel— insulated copper wire wound on a nylon form. Both coils 
are 1.5 mm long, but one has 120 turns with 1.5 mm i.d. and 

3 mm o . d . , while the other has 130 turns with 2 mm i.d. and 

4 mm o.d. Their twisted leads are shielded with a copper braid 
sleeve and housed in a 4 mm o.d. x 10 cm long Pyrex tube holders. 
To allow the coils to probe the upstream side of the anode, the 
holders have right-angle bends 2 cm from the coil. The coils 
are potted in epoxy to protect them from the discharge and the 
tube ends are sealed with epoxy to preclude out-gassing. The 
shielded leads from each probe are passed through Tygon tubing 
to a feed-through at a side port of the vacuum tank and the 

end of the Tygon tube is vacuum sealed to the Pyrex holder. 

The coils are terminated with 11 msec passive integrators^^ 
at the input to Tektronix 555 oscilloscopes with Type 1A7A 
plug-in pre-amplifiers, which are used to monitor the output. 

To calibrate the probes, a known magnetic induction field 
is created in the discharge chamber by replacing the thruster 
cathode with an aluminum post, which is connected to the anode 
by an aluminum plate. From the known geometry of this short 
circuit, the magnetic field is calculated for various discharge 
currents. When the probe output voltage is measured over a 
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range of currents and probe locations, the calibration con- 
stant is established . The calibration constants for the 120 

and 130 turn probes are 1.89 (10 ) and 3.41 (10 Webers/ 

2 

(m -volt) , respectively. The uncertainty in the calibration, 
determined from scatter in the data, is less than 5%. This 
calibration is verified in the actual MPD discharge. The coil 
probes are tucked behind the anode near the chamber backplate, 
beside the insulating ring, so that all current flows down- 
stream of the probe during a quasi-steady discharge. A fur- 
ther check on the calibration, made by locating each probe 
separately at the same place in the discharge chamber, shows 
equal magnetic field for the same current . The risetime of 
the signals from the magnetic probes are compared to the rise- 

(47) 

time of the Rogowski loop that surrounds the accelerator. 

The response of 10 psec which is measured assures that the 
probe signals accurately represent the changing magnetic field 
during the discharge. 

A triple overlay of signals from the 120 turn probe, in 
a 16 kA total current discharge with 24 g/sec injected argon 
mass flow, is shown in Fig. 4.1. The coil is located 2 mm 
behind the anode, 1 mm radially inward from the lip. The 
signal rises at initiation of the discharge and is quasi- 
steady after less than 200 psec. The quasi-steady period 
lasts about 650 psec. At approximately 900 psec, the dis- 
charge extinguishes. A baseline overshoot follows the dis- 
charge as a result of "passive integrator droop". In the 
figure, the integrated output is about 1300 pvolts at 400 psec, 
the time when magnetic field data are usually recorded. This 
output corresponds to a magnetic induction field of about 
4.5(10 ) Webers/m . The discharges are probed throughout a 

region from 2 cm from the anode, to within .2 mm from the anode 
surface, on a 1.5 mm by 1.5 mm grid. 



FIGURE 4-1 
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For the essentially axisymmetric discharges studied, 

the magnitude of enclosed current at each probe location is 

, , (57) 

derived from magnetic field data using Ampere's law. 

Since the grid spacing is relatively fine, the results are 
linearly interpolated and enclosed current contours, i.e 
lines of equal current downstream, are drawn. An example for 
a 16 kA current with 24 g/sec mass flow is shown in Fig. 

4.2. The contours, which also represent current streamlines, 
indicate that current bows downstream, and only 2 kA of the 
16 kA total current enters the anode on the downstream face. 
The dashed lines are extrapolations of the measured current 
contours to the anode surface (the coil size would not per- 
mit probing more closely) . The extrapolation may be mis- 
leading in the region behind the anode, near the chamber in- 
sulating ring, because the Hall parameter (the ratio of elec- 
tron gyro frequency to Coulomb momentum- exchange collision 
frequency) v ' based on measurements described below, is 
greater than unity. The electrons may therefore turn down- 
stream on their way to the anode. 

The local anode current density is calculated from the 
enclosed current contours by dividing the current between 
contours with the anode area that absorbs it. The distribu- 
tion of local anode current density for the 16 kA arc current 
and 24 g/sec argon mass flow is shown in Fig. 4.3. The ab- 
scissa represents a linearly rolled-out anode with midplane 

at 0.0 cm. Like the heat flux distribution, the current dens- 

2 

ity peaks at the lip. It has a maximum there of 300 A/cm , 

2 2 
dropping to 80 A/cm in the discharge chamber and to 10 A/cm 

on the anode face. A check on the integral of the derived 
current densities over the anode agrees with total arc cur- 
rent to within 10%. 

When the local anode Current density, j , is scaled with 

a 

J/A the arc current divided by the anode conduction area 
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9 A 

(y\^260 cm^) viz. j = j /(J/A), the change of current dens- 

a a 

ity distribution with arc current and mass flow can be seen. 

When the mass flow is held at 24 g/sec, the distribution 

(Fig. 4.4) remains almost unchanged, for arc currents from 

8 to 22 kA {J 2 ./m from 3 to 20 kA 2 -sec/g) . On the other 

hand, when mass flow is changed from 6 to 48 g/sec for 16 kA 

fixed current {Fig. 4.5), the portion of total current con- 

ducted over the anode face increases: For 20 kA -sec/g there 

is a threefold increase and for J /m 40 kA -sec/g, as much 

as half the current is being carried over the anode face. 

2 2 

However, if the ratio J /m * 10 kA -sec/g, approximately the 
same portion of current is carried at each anode location. 
This is shown in Fig. 4.6, where the current distribution is 
represented over the arc power range from 0.55 to 3.4 MW. 

For all the operating conditions studied, which fall in the 
current and mass flow range just described, the local anode 
current density is peaked at the lip. 


4.3 Electron Temperature in the Anode Region 


Electron temperature, T , is needed to evaluate electron 

(59) e 

enthalpy ' ' 


h = | £«- 
£ e 


(4.2) 


as well as for the reduction of floating potential measure- 
ments to plasma potential. Both the enthalpy and the plasma 
potential are required to evaluate the anode heat flux. 

The electron temperature is measured with double electro- 
static probes. The probe electrodes are 0.7 cm long x 0.007 cm 
diam tungsten wires, mounted 1 mm apart in ceramic capillaries. 
The rest of the probe construction is similar to the magnetic 
probes described above. 
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To assign a temperature to the electrons, their dis- 
tribution is assumed to be Maxwellian . Therefore, from 
the double probe characteristic , 



A 

where I is the saturation current to the probe extrapolated 
to the current axis at zero potential, and (&I/&>V) q is the 
slope at the origin of the increasing portion of the charac- 
teristic.^^ Current flows between the electrodes when they 
are differentially biased by an electric potential. This probe 
current is monitored with a Tektronix P6042 current probe using 
a 555 oscilloscope with Type L plug-in pre-amplifier. Since 
the probe currents are very sensitive to the surface conditions 
of the electrodes, the electrodes are impact cleaned by ion 
bombardment after each discharge. This cleaning is performed 
by pre-filling the tank with argon gas to a pressure of about 
100 mtorr and then biasing the probe electrodes negatively 
with respect to the anode. The uncertainty of the electron 
temperature measurement is estimated to be about 20%, .since 
the reduction of the data suffers from small superimposed 
fluctuations on the probe current, which may be caused by ion 
density fluctuations in the plasma. 

Electron temperatures measured at the anode lip are shown 
in Fig. 4.7. When argon mass flow m is 24 g/sec, the elec- 
tron temperature increases from 0.7 to 2.0 eV, as arc current 
J is raised from 8 to 22 kA. There is a much weaker depen- 
dance on mass flow from 6 to 48 g/sec, when arc current is a 
constant 16 kA. However, for J /ra ==■ 10 kA^ — sec/g, electron 
temperature is approximately constant at 1.2 eV, over nearly 
an order of magnitude of arc power from 0.55 to 3.4 MW. These 
observed electron temperatures can be represented by the 
following relations: 


UT = O^e^p/o.OUj) 

e. r \ j 


m - 2.4- o/ sec. 


(4.4a) 
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^ = I 4- expt 86©**) /] s T — U.W& (4.4b) 

W-T S' I.? -j 3 /^ S’ ok^-sec/a (4.4c) 

"o -e, j —* 

The results show that electron temperature, and hence electron 
enthalpy, depend on operating conditions, except for constant 
J^/m Is 10 kA^ -sec/ g . 

For the arc currents and argon mass flows examined, the 
electron temperature is relatively uniform over the anode 
surface. The electron temperature 0.2 cm upstream of the 
face of the anode and 0.2 cm from the chamber insulator are 
approximately 20% lower than at the lip; and at a corres- 
ponding position on the downstream face it is about 10% 
lower than at the lip. 

4.4 Anode Region Potential 

Potentials in the vicinity of the anode are determined 
from Langmuir probe measurements of floating potential. The 
probe used has a single cylindrical tungsten electrode 0.1 cm 
long and 0.007 cm diam, radially oriented in the discharge 
chamber. Its construction is similar to the double probe 
used to measure electron temperature. To ensure that no net 
current is drawn, the Langmuir probe is terminated with a 
Tektronix P6013 voltage probe with 100 Mft input impedance. 
Signals are recorded on Tektronix 555 oscilloscopes with 
Type G plug-in pre-amplifiers. An example of the floating 
potential in the anode midplane, 1 mm from the anode surface, 
for 16 kA x 24 g/sec arc operation, is illustrated in Fig. 

4.8. Since the anode potential is ground, the probe re- 
sponse is negative. The signal, like the arc operating volt- 
age included in the figure for comparison, rises in about 
10 psec and becomes quasi-steady in less than 200 psec. 

Errors in the floating potential measurement are small since 
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the data are reproducible, the probe location is well known 

and potential gradients are small except in the anode 

sheath. Accuracy is further assured since the probe 

radius is much larger than the Debye length (which may he 

about 10 m near the anode in the MPD discharges under 

. 21-3 

study, where the electron number density is about 10 m 
and the electron temperature is about 1.2 eV) . 

The probe surveys the anode region over a 1.5 mm x 
1.5 mm grid, to within 1 mm from the surface and lines of 
equal floating potential are mapped. An example is shown 
for the 16 kA x 24 g/sec case in Fig. 4.9, where the dis- 
charge chamber insulator and enclosed current contours are 
drawn for reference. The largest potential gradient appears 
at the anode surface and the floating equipotentials suggest 
a maximum in electric field slightly upstream of the anode 
midplane. The equipotentials in the lip region are almost 
normal to the enclosed current contours. This implies that 
momentum exchange collisions are more frequent than electron 
gyrations about magnetic field lines and therefore, current 
conduction near the lip is scalar. 

If floating potential is examined along current stream- 
lines, as shown in Fig. 4.10 for the 16 kA x 24 g/sec dis- 
charge, and the measurements are extrapolated to the anode, 
the voltage intercept is approximately equal to the measured 
floating potential 1 mm from the surface. Therefore, anode 
adjacent floating potential is defined as the floating poten- 
tial 1 mm from the anode. The distribution of anode adjacent 
floating potential over the anode surface, for 16 kA x 24 g/sec 
operation, is shown in Fig. 4.11. The value in the chamber is 
about 10 volts, dropping to nearly 8 volts in the anode mid- 
plane, and rising again to 13 volts along the downstream face. 

Measured floating potentials, , are not equal to the 

plasma potentials, V , at the probe location. The difference 

ir 

is a shift toward cathode potential because the electron 
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thermal flux to the probe electrode surface is much . higher 
than the total flux of the streaming ions. The calculation 
for the shift is shown in Appendix C. The correction term 



(4.5) 


for a plasma with a streaming ion velocity component u' 
normal to the probe axis, is written in dimensionless form 

_av _ o f I /IS? +J] ■ 

7^?)" I JJ (4.6) 


where rrm is the ion mass. The dimensionless correction 
— & V/ (kT e /e) , shown in Fig. 4.12, is not a strong function 
of electron temperature T , but does vary with ion stream- 
ing velocity. The dashed line in the illustration represents 
the correction when ion streaming velocity is smaller than 
the Bohm speed. (^3,64) 



(4.7) 


The ion streaming velocity component normal to the axis 
of the probe electrode assumed for use in the correction is 
shown in Fig. 4.13. This figure is based on Doppler shift 
measurements of plasma velocity on the discharge centerline 
for a 16 kA x 6 g/sec discharge. At the anode midplane, 
along the axis of the accelerator, the ion streaming velocity 
is about 8700 m/sec downstream. However, the flow angle off- 
axis is 15° with respect to the center of the discharge, with 
a virtual origin at the cathode tip. At the anode, the 
plasma is assumed to follow the contour of the surface, start- 
ing at nearly zero velocity but accelerating as it moves past 
the lip and over the anode face. Therefore, the velocity 
component perpendicular to the radially oriented probe elec- 
trode close to the anode starts at zero in the chamber, rises 
in the lip region and then drops to zero again on the face. 
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Since thrust, the product of mass flow and exhaust 
velocity is a function of current only, (^6,67) a fourfold 
increase in argon propellant mass flow rate, from 6 to 
24 g/sec, may correspond to a four-fold decrease in ion 
streaming velocity, from 8700 to 2200 m/sec, at the same 
16 kA current. Therefore, the dimensionless probe cor- 
rection would change from about -3.5 to about -4.5 and the 
actual probe correction for a 1.2 eV electron temperature 
plasma would change from -4 to -5 volts. This compares to 
the 1 volt error due to electron temperature uncertainty. 
Similarly small changes in the correction might appear if 

the ion streaming velocity in the plasma were to be higher 

( 25 ) 

than 8700 m/sec. 1 Consequently, the velocity component 
normal to the probe axis, shown in Fig. 4.13, is used to 
calculate the floating to plasma potential correction for 
all operating conditions. 

The anode fall voltage is derived from the measured 
anode-adjacent floating potential. Both profiles shown in 
Fig. 4.14 for 16 kA x 24 g/sec operation have the same shape 
along the anode. There is a dip in the lip region and the 
fall voltage is roughly half the floating potential.- If a 
less peaked streaming velocity distribution had been assumed, 
then the depression in fall voltage at the anode midplane 
would be more acute. The error bar is large enough to allow 
the possibility that the anode fall voltage may even be uni- 
form over the anode. 

The anode fall voltages for other operating conditions 
are summarized in Fig. 4.15. For a low current of 8 kA and 
a mass flow of 24 g/sec (Fig. 4.15a), the fall voltage is 
nearly constant over the anode surface at 5 volts, except 
near the lip where it drops to 3 volts. As current is in- 
creased for the same mass flow, the voltage changes only 
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slightly, until the current is raised from 16 to 22 kA 
when the fall voltage over the anode face rises a factor 
of two, to 10 volts. 

For constant current operation at 16 kA, Fig. 4.15b 
shows little difference in anode fall as mass flow changes 
from 48 to 24 g/sec . The anode fall voltage is about 5 
volts in the chamber and on the face, but is depressed to 
about 3 volts on the lip. If argon mass flow is halved 
from 24 to 12 g/sec, the voltage doubles, while the shape 
of the distribution is maintained. Halving the mass flow 
again, to 6 g/sec, causes a maximum anode fall of about 
22 volts in the vicinity of the anode midplane. 

If J /m is held constant at about 10 kA — sec/g, over 

an order of magnitude arc power range from 0.55 to 3.4 MW 

(Fig. 4.15c), much less dramatic changes occur. For 0.55 

MW operation (8 kA x 6 g/sec), the profile is flat at about 

6 volts over the entire anode. Then, as arc power increases, 

a voltage dip in the anode midplane begins to develop. By 

3.4 MW (22 kA x 48 g/sec), the anode fall voltage is 5 volts 

in the chamber, 4 volts at the lip and 10 volts on the anode 

2-2 

face. Thus, if J /m 10 kA -sec/g, the anode fall voltage 
changes slightly in the chamber and rises less than a factor 
of two on the anode face, as the arc power increases. 

Comparison of the anode fall voltage distribution with 
the current density over the anode (Fig. 4.4 to 4.6) shows 
anode fall minima at the highest current densities. 

4.5 Convective and Radiative Heat Flux 

The convective and radiative contributions to the anode 
heat flux, which appear in Eq. (4.1), are negligible com- 
pared to the directly measured anode heat flux. Estimates 
that use conservatively assumed values of the properties re- 
quired are made for the 16 kA x 24 g/sec arc discharge. 
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The order of magnitude of the convective heat flux to 

the anode can be estimated if the anode is regarded as a 

nozzle wall. Start from a dimensionless boundary layer 
(68,69) 

equation 


v _ O. © o 

Nvjk = O.oaXo Re R- 


(4.8) 


where Nu, R g and Pr are the Nusselt, Reynolds and Prandtl 
numbers respectively. This relation may be written 


h 3 D_ e 


=, o.oau 




lO-A- 


(4.9) 


where h is the convective heat transfer coefficient 

3 -2 

D is the scale length = 10 m, the anode width 

is the maximum assumed gas velocity 8700 m/sec 

, o ( 70 > 

kp is the plasma thermal conductivity ^1.0 W/m- K 
jjl is the viscosity of argon ^ 3(10 ^ ) kg/m-sec 


is the plasma density ^ n.m. 

^ i r 


P ^ 


is the specific heat of argon at constant pressure 


= 2000 j/kg-°K ( 71 -) 


In the above, n^ and m. are the argon ion number density and 
mass. Near the anode lip the number density is estimated as 
n^ = 10^m ^ , so that 


So = 


rv* m 

v 


- 7(10 s ) 


(4.10) 


The fluid properties are assumed to have thermodynamic values 
corre sponding to an argon plasma in equilibrium at the elec- 
tron temperature of 1.2 eV (14,000 °K) . Therefore, 


R e - lh±_§ P ? zoo 


(4.11a) 
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and 


Pr- _ = O. O<o 

Mp 


N - = 0.0 s(o(aoo)°' %. otf 4 ^ o. c= 


(4.11b) 


(4.11c) 


Thus, 

S' & <-&©"*), (4.12) 

Therefore, the estimate of convective heat flux to the anode 
lip is 

V*"* - ^ (4 . 13) 


This is two orders-of-magnitude smaller than the directly 
measured heat flux to the lip and can be neglected. As a corol- 
lary, convective cooling due to pre- and post-discharge mass flow 
is also negligible. 

{ 12 \ 

Comparison of plasma radiation mechanisms 1 shows that the 
most intense radiative heating results from '‘free-free” Brems- 
strahlung from the 1 cm long x 0.5 cm diameter plasma blob 
situated near the cathode tip. Using the Spitzer formula 



It-CieT 4 **) 


2. 

rv 

V. 




(4.14a) 


For a singly ionized (Z-l) cathode plasma with an ion density 
23-3 

of n^ = 10 m and an assumed electron temperature of 2.4 eV 
{T e = 28,000 °K) , a factor of two larger than measured in the 
flow, the radiative emission is 

e = aoc> yy/om 3 


(4.14b) 
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Therefore, the radiation from the cylindrical plasma volume 
V ol at the cathode tip is 


\/ 


M & j^ s 4 j>o v'jsWfe 


(4.15) 


If the radiation is assumed uniform over the solid angle 
around the cathode tip, the radiative flux density at the 
anode orifice (r = 5 cm) is 




(4.16) 


which is four orders-of-magnitude less than the maximum mea- 
sured heat flux at the lip and two less than the lowest heat 
flux measured (in the discharge chamber) . 

Although the temperature rise of the inside surface of 
the shell anode has been measured, calculation of the heat the 
anode radiates into the plasma requires the temperature rise of 
the anode surface exposed to plasma. The- plasma-exposed surface 
temperature is estimated by assuming that the anode is a slab 
0<y<JL of zero initial temperature ©• - that there is a con- 
tinuous uniform heat flux q over the surface at y - and there 
is no heat flow over the surface y = 0. This constitutes a stan- 


dard heat flux problem, 
quantity 


(73) 


The solution uses the dimensionless 


R. - 


Je>± © 



JL a 


(4.17) 


graphed in Fig. 4.16 as a function of the ratio y/j£. , for dif- 
ferent values of the parameter «3£T/JL 2 , where T is time. 

For the aluminum anode, = 0.1 cm is the shell thickness 
and the heat flux duration is T = 1 msec. The thermal con- 
ductivity is k = 1.76 W/cm-°C and the thermal diffusivity is 
= 0.823 cmysec. Therefore, 
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JeJ/j*. =0.0© U.is) 

and the curve in Fig. 4.1b, for J^T/J£_ " ~ 0.08 is selected. 

At y/JL.- R=0.24. Since the heat flux q at the anode 

2 a 

lip is about 6000 W/cm , the surface temperature on the heated 
side of the slab at the end of the discharge period is in- 
creased by 





uo°c 


(4.19) 


above ambient room temperature of about 20 °C . Thus, the 
maximum surface temperature of the anode at the end of the 
discharge is about 130 °C (which is much lower than the 660 °C 
melting point of aluminum*' 7 ^). Therefore, the maximum radiated 
heat flux from the anode surface at a temperature Q is 



Cr Q> ¥ o. a w/cm 1 


(4.20) 


TO O 

where CT= 5.7(10 ) W/cm -°K is the Stef an-Boltzmann radia- 

tion constant. Therefore, radiation heat flux to and from the 
anode is negligible. 


4.6 Anode Power in the MPD Arc 

The local anode heat flux and power are calculated below 
from measured plasma properties. The calculation uses a model 
of the anode heat flux based on electron transport within the 
MPD arc. 

As electrons progress through the MPD discharge, they 
acquire energy from the electric field. A portion of this 
energy is surrendered in collisions with the ions and neutrals 
m the plasma. Just before being absorbed by the anode, the 
electrons suffer a last energy exchange collision, after which 
any energy they gain from the electric field cannot be dissipated 
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in the discharge. This energy is added to the electron 
enthalpy prevailing at the point of the last collision and 
is delivered to the anode as heat. 

In the conventional high pressure or atmospheric arc 
discharges, where Eq. (4.1) applies, the distance between 
the position of the last energy exchange collision and the 
anode sheath is small enough and the prevailing electric 
field weak enough that the field-acquired energy electrons 
deliver to the anode closely corresponds to the anode fall 
voltage. However, in the MPD arc the distance between the 
anode sheath and the last energy exchange collision is suf- 
ficiently large, and the .local electric field large enough 
that the anode fall is substantially augmented by a portion of 
the interelectrode field (the field between anode and cathode 
sheaths) . 

The last electron energy exchange collision before strik- 
ing the anode is presumed to occur when the distance measured 
along current streamlines between energy exchange collisions 
X* equals the distance s along a current streamline to the 
anode surface:'*' viz. 

X = S ( 4 . 21 ) 

This energy exchange displacement, x , is derived for the MPD 
arc in Appendix B-l. The displacement is a linear function 
of the electric field and the energy exchange mean free time: 

It increases if potential gradient increases and increases as 
the energy exchange cross section decreases. The influence of 
magnetic field appears in the Hall parameter ^ 58 ^ : A large 

I It must be acknowledged that this is an approximation to the 
actual energy exchange displacement which would vary along 
current streamlines depending on plasma properties. Therefore 
7^ should ideally be calculated as an integral of differential 
energy exchange displacements along the current streamline. 

The data, however, does not warrant the more accurate computa- 
tion. 



79 


Hall parameter indicates increased electron residence on the 
azimuthal magnetic field lines and therefore a decreasing net 
displacement. The potential difference V between s =A_ and 
the anode surface represents field energy delivered to the 
anode as heat. Thus, the anode heat flux model becomes 

(4 - 22) 

where convective and radiative contributions have been omitted. 
In the MPD arc, the energy exchange process is dominated by 
elastic energy exchange collisions between electrons and ions 
(Appendix B-II), The evaluation of the energy exchange dis- 
placement requires electron number densities as well as elec- 
tron temperature and electric and magnetic fields determined 
from probing the plasma. 

The radial electron density profiles used were determined 

previously from a reduction of the spatially unresolved profiles 

( 75 ) 

of the Stark -broadened H^and H^lines. ; The profiles, shown in 

Fig. 4.17 for the 16 kA x 24 g/sec discharge, peak at about 
22 - 3 

10 m along the discharge chamber axis, but drop an order of 
magnitude at the anode orifice radius (5 cm) . 

The measured plasma properties are used to determine the 
dimensionless displacement 



(4.23) 


shown in Fig. 4.18 for the 16 kA x 24 g/sec discharge, as a 
function of the dimensionless streamline co-ordinate 


* 


S E 



(4.24) 


where L is the streamline length from the anode to the cath- 
ode. The error bar in Fig. 4.18 represents uncertainties in 
the electron density by 100%, electron temperature by 20%, and 
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electric field by 20%: The calculation is very strongly 

dependent on the local plasma properties. The results for 
the five current contours listed in the figure fall on about 
the same curve 

* „ 

X !.0,05(s ) (4.25) 

The last energy exchange collision is presumed to occur at 
"X =1 or X=s,i.e. where s =£0.17. Along the L=6.8 cm long 
8 kA current streamline for the 16 kA x 24 g/sec discharge, 
s*-'0.17 represents a distance 

5=Ls*= t,.8(o,q')'\.£cm ( 4 -26) 

from the anode, about one-fourth the anode orifice radius. 

The locus of points where X* = 1 is the "anode power 
surface". It is shown by the circles in Fig. 4.19, where, 
for reference, the discharge chamber, anode, cathode, en- 
closed current contours and floating potentials are drawn. 

The voltage contribution to the anode heat flux, v_ , is 

o 

the potential difference between the anode power surface 
and the anode. The contribution is assessed from the float- 
ing potentials along current streamlines, shown in Fig. 4.20 
for the 16 kA x 24 g/sec discharge. These floating poten- 
tials are corrected to the plasma potentials illustrated in 
Fig. 4.21. Figure 4.21 also includes a plot of the relation 
"X* = 0.05 (s*) 5//3 . For the 16 kA x 24 g/sec discharge, the 
voltage contribution to the anode heat flux, along the 8 kA 
streamline, is V g = 11 volts. Thus, the 4 volt anode fall 
along the same streamline, is augmented by about 7 volts 
of interelectrode potential. 

The anode power is calculated by integrating the results 
of the anode heat flux model, Eq. (4.22), over the anode sur- 
face. The calculated anode power fraction for the 16 kA x 
24 g/sec discharge is 20% of the 1.3 MW arc power. This 
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agrees with the directly measured anode power fraction for 
the same operating condition. If the contribution from the 
interelectrode potential had been neglected, only half the 
anode power would have been accounted for. 

The energy exchange displacement analysis can be applied 
to operation at other arc powers with the same j 2 /m£lO kA 2 - 
sec/g.' The measured electron temperatures, current densities 
and anode fall voltages are used. The interelectrode potential 
is scaled along current streamlines for the 16 kA x 24 g/sec 
case and particle number densities are scaled linearly with 
argon mass flow (assuming equal degrees of ionization, which is 
reasonable since J /m is constant) . The calculated results are 
shown by the circles in Fig. 4.22 where the anode power fraction 
is built up as the sum of the contributions* from electron tem- 
perature, work function, anode fall and voltage from inter- 
electrode potential o From the figure, the voltage contributions 
represent more than half the anode power fraction. Figure 4.22 
shows the agreement between the indirect estimates, using plasma 
properties, and the direct measurements of anode power (solid 
line) over nearly an order-of-magnitude range of arc power from 
0.55 to 3.4 MW. 

During multi— megawatt operation, when the number densities 
in the MPD discharge are high, the anode power surface is close 
to the anode and the commonly employed theory based on anode fall 
and the model presented here based on energy exchange displace- 
ment give nearly identical results. This agreement is a conse- 
quence of the decreased contribution of interelectrode potential 
associated with potential gradients in the discharge, which are 
small compared to the gradients in the anode sheath. The re- 
sult suggests that for a given current and mass flow, arc 
geometry may be adjusted to reduce the anode power fraction to 
a minimum determined by the anode fall voltage required to 
carry the current. 

1 * ' 

The work function for aluminum is about 4 eV. This is 
the same as for almost all metallic conductors. ( 7 6) 
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Explanation of the decrease of anode power fraction with 

2 * 4 

increasing arc power, at the same value of J^/m, now follows 
on the same basis. The electron temperature and work function 
are constant. Therefore, their contributions to the anode 
power fraction decrease with increasing arc voltages at 
elevated arc powers. In contrast, the voltage contributions 
from interelectrode potential and anode fall, whose sum is the 
largest portion of the anode power fraction, change with arc 
power. However, the change of their contribution can be 
assessed because the current distribution is constant if J /m 
remains at 10 kA -sec/g. The anode fall voltage portion of 
the anode power fraction decreases by half, from 0.08 at the 
low arc power of 0.55 MW to 0.04 at 3.4 MW: whereas the inter- 
electrode potential portion, which dominates for the relative- 
ly low power arcs, decreases twice as fast from 0.15 to 0.04, 
as densities near the anode increase. The result is a quickly 
changing interelectrode potential contribution to anode power 
fraction, which falls from its dominating role for low arc 
powers to much less prominence during high-power operation. 

The rapid decrease of the interelectrode potential contribu- 
tion, due to increased number densities in the MPD arc, causes 
the decrease of anode power fraction as arc power increases. 

This model also deals satisfactorily with the variation 
of anode heat flux with mass flow at fixed arc current. An 
example is shown in Fig. 4.23 for 16 kA discharges with mass 
flows above 12 g/ sec (J /m ■< 20 kA -sec/g), where ablation 
and mass recirculation are negligible. In this range, anode 
fall voltage, electron temperature and work function are 
virtually independent of mass flow, while the interelectrode 
potential contribution to anode power decreases with mass 
flow in the manner calculated for the constant J 2 /m case above. 

The total anode power fraction thus decreases with mass flow, 

2 * 

hence increases with J /m, as plotted. The interelectrode 
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field contribution which represents over 25% of the anode 
2 2 

power at J /m ^ 10 kA -sec/g, can be neglected for lower 
2 

values of J /m. The close agreement, shown in Fig. 4.23, 
between the anode power results based on thermocouple data 
and those calculated from plasma properties supports the 
validity of the energy exchange model of anode heat flux 
in the MPD arc. 
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CHAPTER 5 

CONCLUDING DISCUSSION v 

2 

5.1 Anode Power and J /m 

The measurements of anode power in the MPD arc (Chap- 
ter 3) showed a decrease of anode power fraction, for all 
2 

values of J /m, as arc power increases; and this decrease 
was explained (Chapter 4) , for arc operation with 

2 . r*, 2 

J /m =• 10 kA -.sec/g. The explanation asserts that the de- 
crease is due to a diminished interelectrode potential con- 
tribution to the energy electrons deliver to the anode as 
heat, a consequence of increased particle number densities 
near the anode associated with increased mass flow rates 
during high power operation. Now, if the anode power is ex- 
pressed as the ratio of anode power to arc current P /j as 

jHl 

has been done in evaluating steady arc performance, (11,32) 

2 . 

and mass flow is expressed by the ratio J /m which is valid 
in the thruster studied if J /m < 20 kA -sec/g; the result 
in Fig. 5.1 is derived. The figure shows a spread in the 
data that lies outside the representative error bars shown 
for 16 kA total current, but the results do cluster around 
a single curve. This aggregation implies that changes in the 
proportion of the energy electrons carry to the anode as heat 
depend much more strongly on mass flow than on discharge cur- 
rent. Furthermore, P,/J increases when J /m is raised toward 
20 kA -sec/g, i.e. when mass flow is decreased, suggesting 
that for all operating conditions examined, the explanation 
of anode power in the MPD arc is consistent with the data. 
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5.2 Anode Design 

The measurements and the anode heat flux model that 
have been discussed can be used to suggest design changes 
for the MPD arc which may reduce anode power fraction over 
the entire range of arc operation. However, it must be 
acknowledged that any such changes in design may drastical- 
ly alter other parts of the complex arc process. 

Local anode heat flux may be reduced by increasing num- 
ber densities near the anode which decreases the inter- 
electrode potential contribution to the anode power. This 
contribution may ultimately be reduced to zero, if the anode 
adjacent number density is raised to that in the high pres- 
sure arcs ^ ^ or greater. One possibility to increase the 

particle densities would be to supply a portion of the propel- 
lant mass through a porous anode. The increased heavy particle 
density adjacent to the anode should then permit a larger por- 
tion of the energy electrons gain in the interelectrode field 
to be absorbed by the heavy particles and used in the dis- 
charge. In this way otherwise "wasted” power adjacent to the 
anode may be usefully employed for the production of thrust. 
Consequently, unless the arc operating voltage increases more 
quickly than the interelectrode potential contribution is re- 
duced, the anode power fraction will decrease. 

Measurements in the discharge showed that the anode fall 
voltage is generally lowest at the lip, where current density 
is a maximum. The current conduction area of the anode was 
therefore restricted to a 1 mm wide band at the lip just up- 
stream of the midplane. For 16 kA x 24 g/sec arc operation, 
the terminal voltage dropped by 25% while the anode adjacent 
floating potential fell by 10% compared to that of a "stan- 
dard" anode. These results suggest that less power is being 
consumed by the anode. They also confirm the preliminary 
measurements of anode heat flux, on a thruster with current 
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conduction at the anode over a 2 mm wide band at the lip, 

which demonstrated anode power fractions generally half those 

for a thruster with a fully exposed anode, at the same opera- 

(99) 

ting conditions. On the basis of these measurements, it is 
suggested that current conduction at the anode be limited to 
a small band at the anode lip. 

For a thruster with a given geometry, the choice of pro- 
pellant might also reduce the anode power fraction, if it has 
a large electron-heavy particle energy exchange collision cross 
section compared to argon. The anode power surface may then 

move closer to the anode, thereby decreasing the anode power 

(77 ) 

fraction. A suggested propellant is lithium, ' ' although 

this major change may presage the onset of other modes of 
power consumption and operating difficulties. Lithium has 

cross sections for energy exchange at least four-orders-of mag- 

( 78 ) 

nitude larger than argon due to its low-lying excited states. 
However, for the same reason, it is also notorious for large 
frozen flow losses, which may exceed any benefits to be gained. 
Nevertheless, although it is highly reactive, lithium is easy 
to store as a solid and can be liquified as it regeneratively 
cools the anode to regain some of the residual anode power. 

For those reasons, research on a multi-megawatt, lithium-fed, 

MPD thruster would not only be instructive but may prove bene- 
ficial to thruster improvement. 

Reduction of the anode power by decreasing electron tem- 
perature and anode material work function is more difficult. 

The electron temperature is established by electron collision 
processes in the discharge and it tends to hover around 1 to 
2 eV. Although propellants with higher inelastic energy ex- 
change probabilities may be used to replace argon to decrease 
the electron temperature, the frozen flow losses may outstip 
any efficiency gains from anode power reduction. To decrease 
the work function contribution, different anode materials are 
required. Presently, the work functions for good conductors 
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all have values of about 4 eV t Improvements in them 

require advances in materials engineering. 

5.3 Conclusions 

This investigation has focused on the power consumed 
by the anode of the MPD arc. A series of direct thermo- 
couple measurements of anode heat flux verify that anode 
power fraction decreases with increasing arc power, as 
suggested by previous indirect experiments. 

Complementing the thermocouple measurements are probe 
studies of the plasma properties in the vicinity of the anode 
Their results are employed in an anode heat flux model modi- 
fied for application to the MPD arc. In the proposed heat 
flux model, a key ingredient is the potential difference be- 
tween the anode and the location of the last effective energy 
exchange collision between electrons and ions before the elec- 
trons are absorbed by the anode. For the experimental con- 
ditions examined in this study, this last energy exchange 
collision is located outside the anode sheatho Therefore, 
the energy which is added to the electron enthalpy and work 
function terms involves not only the anode fall, but an ad- 
ditional contribution from the potential fall outside the 
anode sheath. The difference between the commonly employed 
anode heat flux model and the one suggested here is simply 
this additional voltage contribution. 

The anode power of the MPD arc, calculated on the basis of 
the energy-exchange model, agrees with the results of the 
direct measurements. It explains the decreasing fraction of 
input power consumed by the anode, as arc power is raised, on 
the basis of increased propellant mass flow, since the theo- 
retical location of the last energy-exchange collision approaches 
the anode faster than the local electric field outside the anode 
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sheath increases. Ultimately, the densities may be large 
enough that the relevant voltage contribution can be approxi- 
mated by the anode fall, as for conventional electric arcs. 

The measured decrease of anode power fraction with in- 
creasing arc power suggests that raising arc powers to multi- 
megawatt levels makes a larger portion of the input power 
available to the discharge for the production of thrust, and 
therefore raises thrust efficiency of the accelerator. 
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APPENDIX A 

SHELL ANODE HEAT CONDUCTION 


It is desired to relate the measured temperature rise 
on the inside surface of the shell anode due to a quasi- 
steady discharge pulse to the average heat flux that pro- 
duces it. The problem is solved under the assumptions 
that the heat flux is constant throughout the discharge 
period and is zero at other times, and that the anode 
geometry is locally flat. 

Suppose that the anode can be examined locally as a 
slab of thickness and that the initial temperature is 
uniformly zero. The applicable heat conduction equation 
in one dimension is 

= jfe. m 

3 3-t (a.i) 

where , c, and k a are the density, specific heat, and 
thermal conductivity of the anode material; t is time and 
y is distance through the slab with one face at y = o and 
the other at y = & ; and ^ is temperature. 

The problem is linear and its solution is a two part 
sum. First the heat flux problem for a constant uniform 
heat flux q supplied continuously from t = o to the face 
at y = JL is solved. Then the solution is used for a con- 
stant heat flux - q of equal magnitude but opposite sign, 
which starts at t = T later. When these two are summed, 
the result applies to a rectangular heat flux q which 
lasts from t = o to t = T: 
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The geometry of the problem is 



and the initial and boundary conditions are: 
initial uniform zero temperature 




no heat flux across the boundary at y = o 

it ^(o„iVo 

■t 

heat flux over the boundary at y = £ is a single, 
square pulse 


J? 30 

\dy 


r 

U,t) = ^ hc^ - w(^_-n • 
* 


where H is the Heaviside function. 

The uniform supply of heat to one side of the slab 
at y - with no loss of heat through the other side at 
y - o is a standard heat conduction problem and its solu- 
tion appears in Carslaw and Jaeger. 
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where the thermal dif fusivity is 


^ * 


When the heat -q is supplied at t = T, 

, t >T , 

Ik* m y 

Now, since the problem is linear, the temperature 
©Cy) at times t>T is found by adding Qi^y) and ©2^ ' 


thus 


qT ol[z y&£ L L - 

~gi L ^ ^ 


e >?(- 


1 

L ) c ^[ n -fi) , 1>T 
J J - 


and the temperature on the unheated side y = o can be 


written 


i-T 




(A. 9) 


t>T 


The first term represents the final temperature of the 
slab while the second is the transient heating correc- 
tion. 
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I f the relation is written as 


©00 = 


t>r 


(A. 10) 


then represents the temperature rise per unit 

heat flux and so the heat flux can be calculated from 
the temperature rise as 





(A.ll) 


These computations are performed for the shell 
anode for temperatures measured by thermocouples 
mounted on the inside. The calculations are programmed 
for use on the IBM 360/91 computer at the Princeton 

University Computer Center. The series is carried to 

-7 

terms of relative order 10 of the preceding sum. 

To establish how well the heat flux results de- 
rived from the above local, 1-dimensional formulation 
represent a nonuniform heat flux over the anode, it is 
instructive to compare Eq. (A.l) with the more general 
2-dimensional heat flux equation: 


OC, 3 ©- 

" at j 


9L& , 

cJX*- 


L- 

where x represents a coordinate along the 
shell . 

If this equation is integrated through 
thickness 


(A. 12) 

anode 
the shell 



(A. 13) 
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or, since no beat flows at y = o. 


*$ c at = 


3 X © 

3 ^ 




This relation can be approximated by 


(A. 14) 



or 



(A. 15a) 


In 



this 


form. 



represents the correction to 


(A. 15b) 


the 


one-dimensional problem. 

To assess the effect of this correction, a finite 
difference technique is used to evaluate the ratio. 


If A. . 

time j 
elapsed 



represents the temperature at location i 
, and the separation between points is AX. 
time is & t, the ratio becomes 


(A. 16) 
and 

, whi le 
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For arc operation at 16 kA x 24 g/sec, the most severe 
temperature gradients appearing on the aluminum anode are 
in the lip region, near the anode face. At the time the 
data are examined, 20 msec after the discharge. 



2. 2PQo~ 5 ) 

2 7 (o.g) (o.2S) a 


2, - 2.Qq) » 20 

1=)- o 


O.Z. 


(A. 18) 


Thus, at the worst location, the heat flux may have an error . 
of 20%. Elsewhere, the error is less than 10%. It is to this 
approximation that the one-dimensional heat flux relation, 

Eq. (A. 11), is used everywhere on the anode. 
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APPENDIX B-I 

THE ENERGY EXCHANGE DISPLACEMENT 


The purpose of this appendix is to calculate the mean 
displacement along current streamlines between energy exchange 
collisions of the current carrying electrons and the other 
(heavy) species in the cylindrical MPD plasma. 

Consider a "sample" electron drifting from the cath- 
ode to the anode in an ax i symmetric electric field 

IE = (B-i) 

and an azimuthal magnetic induction field 




(B.2) 


( 


where e^, e^ and e z are unit vectors in a circular 
cylindrical coordinate system. 

The motion of the electrons is described by 


~ e ( e + x ^) _ "K (\r — (B . 3a) 

a * 

or 


i r 





lB ) ft - r. ) 

a 1 i 


(B. 3b) 


where yr is a position vector,}^ p , is the momentum- 

® 1 

exchange collision frequency between the electron (e) and 

another particle (j ) t and e and m are the electronic c ha roe 

e -* 

and mass respectively. To simplify the formulation, the 
sample electron is replaced by an "average" electron, such that 
j=e is not considered. 
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To simplify Eq. (B.3b), it is assumed that particles 

of each species are in thermal equilibrium with one another. 

(821 

Therefore, they move at a mean thermal speed * 

v - ' (B.4 

tk / Tdm 

where m is the particle mass and k is Boltzmann's constant. 

Since m.^rn , the ratio of the mean thermal speeds of species j 
j ' e 

to the electron mean thermal speed 

~ 1 

even if T. NT ; i.e. ions and neutrals move much more slowly 
3- e . • 

than the less massive electrons. Therefore, \ry?fcTj 

and so 

ir -ir = ( B . 6 

$ 

from which the equation of motion for the average electron 
becomes 

r - -^(E - rxB)- rZ vj <b.7 

* 

The radial, azimuthal and axial components of this equa- 
tion with an electromagnetic field specified by Eqs . (B.l) 

and (B.2) are 


radial 


*1 s S ^ 


-rJ 


azimuthal : 




axial : 


(B.8c) 
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From these equations the radial and axial motion ( i . e . 
toward the electrodes) is independent of any azimuthal 
motion (although the inverse is not so) . They can be 
written as 


4 d r~- 

m e c 


• 2 hT _o 


^ — z- 21 =o 

n-.- a C - £ 


(B . 9a) 

(B.9b) 


If it is further assumed that the ionization is 
above a few percent, the momentum - exchange is dominated 
by long range Coulomb collisions between electrons and 


ions (i) . 


(84) 


Thus , 

Z V> p s v p 

— e.u 


J d 

and Eqs. (B.7)can be written in the form 


(B. 10) 



where 




\ 


^3 




i -*n3- 

_o. 




m e 

_ <-*=< 


€Lw 



(B. 11) 


( B. 12a) 


B. 12b) 


B . 1 2c) 


are the electron mobility, Larmor radius and Hall para- 

. (58,62,83) 

meter. v 

Now, the integral of Eq. (B.9) is to be found over the 
mean free time between energy exchange collisions 

^ E l 




(B. 13) 



where n^ is the number density of the field particles 
which exchange energy with the electron and is t * ie 

collision cross section for the energy exchange mechanism 

r P'hn e; . 


= 



(B . 14a) 



(B c 14b) 


is the change in the electrons position between energy ex- 
change collisions. This result is valid if the electron 
mobility ,Larmor radius, Hall parameter and electric field 
are constant over the distance between energy exchange col- 
lisions, because the motion is then determined independant- 
ly of energy exchange. To within the accuracy of the experi- 
mental measurements in this study, this is the case in. the 
MPD discharges examined here. 

The position change corresponds to a distance 





(B „ 15) 


which follows a current streamline since electrons carry 
the current. Hence, the energy exchange displacement along 
a current streamline is 
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The displacement ^ is a linear function of the elec- 
tric field and the energy exchange mean free time: It in- 

creases if potential gradient increases and increases as the 
energy exchange cross section decreases. The Hall parameter 
represents the influence of magnetic field: A large Hall 

parameter indicates increased electron residence on the 
azimuthal magnetic field lines and therefore a decreasing 
net displacement. 
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APPENDIX B-II 

ENERGY EXCHANGE COLLISIONS IN THE MPD ARC 


The event which dominates the energy exchange in the 
MPD arc is the one with the largest collision frequency 

E 

■>ih = \O eh \ t (B . 17 ) 

where , the heavy particle number density, may represent 

ions and atoms. 

The energy exchange processes occurring between the 
electrons and the heavy particles are neutral and ion excita- 
tion, neutral and ion ionization, and elastic energy exchange. 
The energy exchange process which dominates in the anode region 
of the MPD discharge is assessed using energy averaged cross 
sections representing collision frequencies f nv. These energy- 
averaged cross sections are estimated below- for an argon plasma 
with a 1.2 eV electron temperature: 

Q^ lon : Ionization of neutral argon atoms by electron col- 

lision. 

The threshold for this event is 15.75 eV, the 
ionization potential of argon, which is well above the 
mean electron energy in the plasma. From Appendix D-I , 
the effective inelastic cross section between electrons 
with a Maxwellian distribution whose mean energy is 
below threshold, and heavy particles which ‘are not 
moving is 


< Q r > * 4!, 












(B.18) 


where E q is the threshold for the event and is 

the energy where a linear approximation to the increasing 
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Q 


E, ex 
ei 


Q 


part of the actual cross section is assumed to have 

a value equal to the maximum cross section cr ^ ^ 

eh' 

For the ionization of the neutral argon atom, the 
maximum cross section, 4{1CT 20 ) m 2 , occurs at about 
16eV. (85 > 

Therefore , 




(B, 1 9a) 
(B.19b) 


E, ion 
ei 


E, ex 
ea 


2 Excitation of argon ions by electron collision. 

Excitation events which populate excited states 

of the argon ion from the ground state have thresholds 

about 17 eV . The maximum cross section is on the order 

of 1.2(10 -20 ) m 2 for 21 eV electrons/ 86-89 ’ Therefore, 

the effective cross section in the 1.2 eV electron 

temperature plasma is about Q E ' ex = (10“ 27 ) m 2 . 

e 1 

: Second ionization by electron collision. 

This event has threshold at 24 eV and a maximum 
cross section of 1.2 (kT 20 ) m 2 for 79.4 eV electrons . (85 ^ 
Its effective cross section is about Q^*ion_ 9(10 _28 )m 2 a 

Excitation of neutral argon atoms by electron 
collision . 


Although the energies of threshold and maximum are 
well documented and resonance energy levels are available, 

the absolute magnitudes of the cross section are not well 
known and can only be roughly estimated. The cross 
section assumed here is conservatively large. It rises 
from an 11 eV threshold to a maximum at 22 eV. At the 
peak it is assumed to be 4(10 -22 ) m 2 , two orders-of- 
magnitude smaller than the maximum cross section for 



no 


ionization of the atom, Equation (B. 17) gives 

n E, ex -3 / -j ,->—26 , 2 

Q ea = 3(10 ) m . 

g p 

Q ei : Translational energy lost in elastic, Coulomb 

momentum-exchange collisions between ions and elec- 
trons. 

This event does not populate ion excited states 
but it is and provides a mechanism way for electrons 
to share energy with ions. The fraction of the inci- 
dent energy lost by electrons in these elastic en- 
counters is 2 (m/m i ) of the incident energy . ^ 91 ' 92 ^ 
Therefore, the energy the electrons have is effective- 
ly exchanged in about 35,000 encounters. For the elec- 
tron temperature of 1.2 eV and an electron number dens- 
21 -3 

j.ry of 10 m in the anode region, {see Fig. 4.17), 
the Coulomb momentum-exchange cross section, from 
Eq.. (4.36), is about 1(10 . Therefore, the energy 

exchange cross section is about 4(10 ) m . This is 

larger than the other cross sections for energy ex- 
change, but the magnitude decreases by about 25% with 
each 10 % increase in electron temperature. 

The cross sections above are now used along with the esti- 
mates of the heavy particle number densities to calculate the 
energy exchange collision frequencies. First the heavy particle 
number densities are related to the electron density in the dis— 
oharge . For the quasi— neutral plasma, the ion and electron num- 
ber densities are equal 


Y\ 


n • 


(B.20) 


and if the ionization fraction, the ratio of ion to ion-plus- 
neutral number density, is oc ; the neutral argon number density 
is 


rv 





) 




( B . 2 1 ) 



Ill 


Thus , the collision frequency for electron-ion energy ex- 
change is 


,(*) _ 


n _ Ol • v 


<=kj£- 


whereas for electron-atom collisions it is 


(B. 22) 


- o (B.23) 

where ( ^) is the type of event. For the 16 kA x 24 g/sec 
discharge, whose electron temperature is about 1*2 eV, in an 
estimated ionization fraction 04 - 0,9. Therefore, 
the ratio (1 - «6) /oc =* 0 . 11 . -The collision cross sections 
above, used in the collision frequencies f ^ andp ^ 

61 6 a. 

from Eqs. (B.22) and (B.23), then show that multiple elastic 
Coulomb momentum-exchange collisions Q ^ dominate the 
energy exchange process in the MPD arc. The multiple elas- 
tic collisions would still exceed the other events if the 
ionization fraction were 0.01, two orders-of -magnitude smaller. 
However, if the electron temperature increases to the neigh- 
borhood of 5 eV, the inelastic exchange mechanisms will dominate 
the energy exchange process. 
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REDUCTION OF FLOATING POTENTIAL 
TO PLASMA POTENTIAL 


Although floating potential measurement is straight- 
forward, the magnitude of floating potential does not equal 
the plasma potential at the probe location. The difference 
is a shift toward cathode potential because the electron ar- 
rival rate to the probe electrode surface, due to their ran- 
dom thermal flux, is much higher than that of the ions. The 
shift is calculated by assuming that the probe is immersed 
in a streaming plasma. Since the probe radius is less than 
the Larmor radius,^ magnetic field effects are neglected . ^ 

An electron current toward the probe exists in the plasma 
because electron mobility is much higher than that of the 
ions. The res*.jlt is a negative sheath which repels electrons 
and attracts ions. However, to attract sufficient ion cur- 
rent to maintain the no net current condition, a "pre-sheath” 
is created between the plasma and the actual probe sheath. 

The pre-sheath is quasi-neutral, but has an electric field 
which accelerates ions. In a quiescent plasma with a Maxwellian 

electron distribution, Bohm has shown that the pre-sheath 

( 63 ) 

accelerates ions to the speed 



where nr is the ion mass, is the electron temperature 

and k is Boltzmann's constant. Thus, the ion current to 

^The Larmor radius due to a magnetic field B is 



eB 


t / SkTe.rrvg \ 

^ ~~ Tt * J 


(C.l) 


where m Q is the electron mass and the electron thermal speed 

For a 1.2 eV electron temperature and a 
magnetic field of 10“2 Webers/m, the Larmor radius is about 
4(10~ 4 )m which is greater than the probe radius of about 10 -4 m. 
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the probe is 


I- = c 


(C.3a) 



(C.3b) 


where Z is the ion charge number, n^ is the ion 
number density and A is sheath area. If the probe 
has a characteristic dimension much less than any macro- 
scopic dimension associated with the plasma and much 

(61) 

larger than the Debye length, then A equals the 

probe conduction area. 

The electron velocity toward the probe is 

-\ 


U 


l 


u.T« 


where m is the electron mass. 


(C . 4) 


Thus, the electron cur- 


rent to a probe with a voltage bias 
respect to the plasma is 


&V ® V - V_ 
p f 


with 


I e 




WT e 



A 


(C.5) 


The pre-sheath will adjust its voltage and thickness 
to maintain the zero net probe current at floating potential 


T - T (c.6) 

Since the plasma is quasi-neutral at the "boundary” be- 
tween the sheath and the pre-slieath. 


rv 


n, 


(C. 7) 


and so 





L 


e-CVp-yh 

kT. 



(c.s) 
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Thus, the correction from floating to plasma potential is 
For a singly ionized argon plasma. 


(C . 9a) 


=- 5 . 4 - 


(C.9b) 


However, in a streaming plasma, the ions may have a 

velocity u which exceeds the Bohm speed /mT' and 

a smaller pre-sheath voltage is required to balance ion 

( 95 ) 

and electron currents. ’ 

Suppose the ion streaming current is 


X- = "2L o* n ♦ A s (C.10) 

os o a 

where u ^ is the component of streaming velocity normal 
to the probe and A is an appropriate ion collection 
area of the probe. If ion current is assumed to be the 
sum of Bohm speed and streaming velocity contribution, 
the current balance becomes 

T- + 1 . - I, < C - Lla ) 




A, + ctLu n; A, 


y ZTCm c 


rv 


- (c.iib) 


V«£T, 


Therefore, 


m=v p _v s 


= -y* ju 

e. 


*/ 


A 

2-Ttgy t 

VcT*. 


r 


A/iS. v A 


rr\ 


(C.12) 
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In Eq, (C.12), ion current is expressed as the sum of 
a "Bohm" current X ; and an "ion streaming com- 

ponent” current When o^u^u^, the Bohm current 

dominates the logarithm. Ion streaming becomes important 
when u » i.e. when the ion velocity component 

normal to the probe exceed the Bohm speed. 

The area for streaming ion current is probably larger 
than the projected probe area normal to the ion stream 
since the probe potential bias will bend the ion trajector- 
ies to draw current. However, it is unlikely that ions are 
collected over the entire probe area. Therefore, half the 
probe area is selected as a compromise for streaming ion 
collection. 



Also, since pf appears in the logarithm argument, its 
value may change by a factor of two without significant- 
ly affecting 

The probe correction for ions streaming toward the 
probe, becomes 





(C . 14) 


and for a singly ionized argon plasma with electron tem- 
peratures in the range from 0.6 to 2 eV (T — 7000 to 
23,000 K) and ions streaming at about 10^ m/sec, the cor- 
rection is 


s-'&.'s 


(C.15) 
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APPENDIX D 

FOLDING OF INELASTIC CROSS-SECTIONS 
WITH ELECTRON DISTRIBUTION 


The electron distribution n (E) and the collision 
cross section (e) for an inelastic event ) between 

the electron and a heavy particle (h) are functions of 
energy E. Schematically they might look like the sketch below: 



The electron temperature for the distribution shown repre- 
sents a lower average energy than the threshold for the energy 
exchange process, but there is a portion of the population 
that is energetic enough to accomplish the inelastic event. 

If the actual energy exchange cross section (E) is 

folded into the electron distribution, an "effective ' 1 
collision frequency and a mean cross section for the event 
can be computed. 
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Suppose energy is written in the form 


where m is mass and velocity and energy dis- 

tributions are replaced by distributions in velocity. 
The mean collision frequency ^ or an ^ n_ 

elastic event may then be calculated from 


(D.l) 





OT I 

^©W 



(D.2) 


where 1 indicates “function of" and where 

Wa* 

n e ( v / e ) and n^V/^) are electron and heavy particle dis- 
tributions in velocity space* is the integral of the 


eo 


electron distribution over velocity space and 
are the electron and heavy particle masses, 
the heavy particles are not moving; 


m 


and 


If 


H 


(D.3) 


Eq. (D,2) becomes 



For a Maxwell-Boltzmann electron distribution about tem- 
perature T 

e 



(D.5) 
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the integral becomes, 


n 


Co 




> 


r n hj\o( 


m 4 




J p 


{- 1 

V 2JteTc_ ) 


C& 

Q 1 




lv/ e \ cW e 


$> 


Let Q eh r be approximated by a cross section that zero 
at threshold and rises linearly with energy with a slope 
equal to that at threshold: 


Qt 


/“O 


IcK 




(fr^l 


where represents the maximum cross section for the 

event, is the threshold energy, and E^ is the energy 
where the maximum cross sectioned is assumed to occur. 
Using the definition of Eq. (D.l), the cross section can 
also be represented as 


c? 


} K - 

n v< - 


(D.6) 


This formula is substituted into Eq. (D.6) to form 


«.o 


^ V'v, f 


oo 


( 2-TCteTw) 


3/, 




(- 

v Z.WX. / 
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,e>* 


. It 

V/e -Mo 

V//; - V 4 l 


iNTeld V/ t 


(D.7) 


(D.8) 


(D.9) 


Since the integral is zero for electron velocities from 0 to 

V eQ , the limits of integration over all velocity space are 

replaced by integration from threshold velocity V . Infinitv 

eo 2 

is retained as the upper limit, since the integrand decays 
exponentially as the square of velocity, the independant 
variable. This use of infinity simplifies the computation 
without seriously effecting the result. 
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The integral is then written as 

cO 


C*> 


)KVj-v/^ |V4\ 


'•* 




y l 

^c*> 


(D. 10} 


To integrate the spherical coordinate system shown below 
is used 



v ‘ n (t a ) °°*iX z) (D ' 1La) 

vl* ,n ^j Vvn (t) (D - Ub) 

^ 3 ) (d. lie) 


Thus 


IvJ. £, 


and 
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(D. 12) 


(D. 13) 


(D. 14) 
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After integration and simplification, Eq. (D.13) 
becomes 

n eo <^ } > « je£ nvv o^ y — ( 


( 


4L 

2.Vc^. 

2 




M- 


2.WX*. / 


(D, 1 5) 


The velocity squared V/ can be replaced by J^E/m 1 
from Eq. (D.l). Therefore, Eq. (D.14) becomes 



(D. 16) 


(D.17) 


and since the electron mean thermal speed is 


Therefore, 

< «‘S,> * 



(D . 18) 


(D. 19) 


Thus, the energy dependent inelastic cross-section is 
replaced by a mean value which depends on the threshold 
energy Eq, the cross-section maximum , the slope of 

the cross section at threshold and the electron temperature . 
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ANODE HEAT FLUX AND POWER: 
COMPENDIUM OF RESULTS 
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